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The study of Majorana fermions is of great importance for the implementation of a quantum 
computer. These modes are topologically protected and very stable. It is now well known that 
a p-wave superconducting wire can sustain, in its topological non-trivial phase, Majorana quasi¬ 
particles at its ends. Since this type of superconductor is not found in nature, many methods 
have been devised to implement it. Most of them rely on the spin-orbit interaction. In this paper 
we study the superconducting properties of a two-band system in the presence of antisymmetric 
hybridization. We consider inter-band attractive interactions and also an attractive interaction in 
one of the bands. We show that superconducting fluctuations with p-wave character are induced 
in the non-interacting band due to the combined effects of inter-band coupling and hybridization. 
In the case of a wire, this type of induced superconductivity gives rise to four Majorana modes at 
its ends. The long range correlation between the different charge states of these modes offers new 
possibilities for the implementation of protected q-bits. 


I. INTRODUCTION 

In metallic multi-band systems, electrons arising from 
different atomic orbitals coexist at a common Fermi sur¬ 
face. Electronic states, either in the same or different 
sites may overlap and mix through the crystalline po¬ 
tential giving rise to hybrid bands. Superconductivity 
is strongly affected by hybridization [lJ-Q. In general 
this has a detrimental effect and can even destroy it at 
a superconducting quantum critical point (SQCP) Q- 
This has been verified experimentally [5|-|7|] in multi-band 
superconductors which are driven to the normal state 
by external pressure which increases the overlap of the 
wave functions and consequently their mixing. Theoret¬ 
ically [§j, |9] this has been verified for the case of a sym¬ 
metric, fc-independent hybridization. However, in many 
cases, as will be discussedin this paper, hybridization 
can be antisymmetric (Toi - pHl and we should consider its 
fc-dependence. Here, we study the influence of odd-parity 
mixing on the superconducting properties of a multi¬ 
band superconductor. We show it has non-trivial effects 
in the superconducting properties. We demonstrate that 
antisymmetric hybridization can enhance superconduc¬ 
tivity m in multi-band systems. We also show that it 
can promote a crossover from a weak coupling, Bardeen, 
Cooper and Schrieffer (BCS) 0 type of superconduc¬ 
tivity to a strong coupling one associated with a Bose- 
Einstein condensation (BEC) of pairs @, 0- 

One of the most remarkable properties of a one dimen¬ 
sional p -wave superconductor was discovered by Kitaev 
fl5| ; he has shown that Majorana fermions can exist at 
the ends of this system. The idea of Majorana fermions 


was introduced by Ettore Majorana in 1937 0. After 
more then 70 years it was proposed, there is still no defini¬ 
tive way to detect it [171-121]. This particle has the eccen¬ 
tric property of being its own antiparticle [22| . In con¬ 
densed matter systems Majorana fermions are emergent 
quasi-particles, topologically protected and that satisfy 
a criterion of robustness for use in quantum computers. 
They are promising candidates to acts as <7-bits [231 ]. In 
the high-energy context, there is a current idea that neu¬ 
trinos may be Majorana fermions [24]. The investigation 
of Majorana fermions is also important in the context of 
fundamental physics. 

It becomes clear that obtaining Majorana fermions is 
of great importance. Motived by this, we propose here 
a new mechanism to produce a p-wave one-dimensional 
superconductor, without the necessity of spin-orbit in¬ 
teractions and an external magnetic field. The two-band 
Hamiltonian we study in this paper represents an effec¬ 
tive model to describe a non-interacting wire deposited 
on top of a bulk superconductor. Assuming that the hy¬ 
bridization between the electronic states in the wire and 
on the bulk is antisymmetric, we show that the super¬ 
conductivity induced in the wire has a p -wave charac¬ 
ter. However, since we are dealing with spinful fermions, 
the pairing we obtain corresponds to the mi = 0 of the 
l = 1 p -wave state. This type of pairing does not break 
time reversal symmetry and we find four Majorana modes 
in the chain, two at each end. Although two Majorana 
can combine to form a standard fermion, these composite 
Majoranas exhibit highly non-trivial properties [25| and 
their charge states have long range correlations along the 
wire. 
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II. ORIGIN OF ODD-PARITY HYBRIDIZATION 

In this section, we discuss the origin of an antisym¬ 
metric hybridization following the work in Ref. [26(. 
We assume, as these authors, that the Wannier func¬ 
tions (13, HU of the s, p, d or / electrons of the solid 
have the same parities as the corresponding atomic func¬ 
tions. 

We consider that hybridization is caused by a peri¬ 
odic lattice potential, v(r). This lattice has inversion 
symmetry, i.e., v(—r) = v(r). The matrix elements of 
hybridization are written as, 

/ OO 

dr'ip*(r' - n)v(f')ipv(f' - f 2 ), 

-oo 

where is the wave function of 1(1') orbital. We can 
take r*i = 0 and r 2 = — f. So, 

/ OO 

dr'ip*(r')v(f')ipi>(r' + r), ( 1 ) 

-OO 

where we can write ip in spherical coordinates [ 1 U as, 


has even parity ill' — l is an even number and odd parity 
if V — l is an odd number. Then, every time we mix or¬ 
bitals in neighboring sites with angular momentum l and 
l + 1, we need to consider odd-parity hybridization. The 
antisymmetric relation in real space is V(—r) = — V(r) 
and in momentum space it is given by, V(—k) = — V(k). 
In a one-dimensional lattice, for example, an antisym¬ 
metric hybridization is given by, V(k) = 14 oc « sin ka 
with a the lattice spacing. An additional important con¬ 
straint is that of time reversal symmetry which implies 
for the Hamiltonian studied here that the antisymmetric 
hybridization 14 is a purely imaginary quantity. 

Notice that the case of antisymmetric V(r) is of great 
relevance for condensed matter physics as it includes the 
s-p 1 hybridization, d-p mixing, relevant for the copper ox¬ 
ides, and d-f mixing that encompasses many rare-earth 
systems, the actinides and their compounds. Antisym¬ 
metric mixing is also an essential ingredient to give rise to 
topological insulating phases in multi-band systems 3_() ■ 

III. MODEL 


MV = 4’im(r,9,(p) = R(r)Yr(9,<P), ( 2 ) 

with R(r) the radial solution of Laplace’s equation and 
Y l m (d, <p) the angular solution, known as spherical har¬ 
monics. The indexes l and m are quantum numbers, 
such that l > 0 and m = —l ,..., 0,..., +1. 

We want to investigate the parity of Eq. |I] This is 
possible by doing an inversion of coordinates: r —> —r. 
Doing this in equation (121) . we notice that the parity of 
the wave function depends on the parity of the Y™, the 
spherical harmonics, which depends on l by the following 
expression, 

Yr(7r-e,7r + cP) = (-l) l Y l m (0,cP). (3) 

So, 

M-V = (~ 1 ) l MVi ( 4 ) 

such that, ip*(—r) = (— l)~ l ip*(r). 

Now we are able to determine the parity of Eq.[TJ Per¬ 
forming an inversion of coordinates f —> —r, we can write, 

/ OO 

dr'ip* (r')v(r')ipi> (r' - r), 

-OO 

doing a change of variable f' = — r ", we get 

/ OO 

dr"iPt(-f")v(f")ip l , (- \r" + r ]). 

-OO 

Finally, using Eq. [4] and doing some simple manipula¬ 
tions, we get, 

Vw(—r) = (—l) l '~ l Vn'(r). (5) 

We can conclude that the parity of the hybridization de¬ 
pends on the difference of the angular momenta of the 
electronic orbitals, I'—l. Since, l is always positive, Vw (r) 


As we mentioned earlier, we focus our attention on a 
two-band system, with an odd-parity hybridization be¬ 
tween these bands, an attractive (inter-band) interac¬ 
tion between the electrons in different bands, and also 
an attractive (intra-band) interaction between electrons 
in only one of the bands. The Hamiltonian of this prob¬ 
lem can be written as, 

H ~ (fk a kcr a k& + e k9k a bka'j 

k,a 

y! {^■abn\ a b^_ k _ a + A ab b-k-a-C1ka^ 

ka 

Y. ( K / ^bbb\.y__ k _ cr + A bb b-k-<jbka^ 

ka 

+ (v k a kcr bk<T + 14 j (6) 

ka 

where a is the spin index that could be “up” (f) or 
“down” ( 4 _), e k ’ b are the energies of the electrons in the a 

and b bands. In an obvious notation, ajjj and b k ^ anni¬ 
hilate (create) electrons in these bands respectively. The 
attractive many-body term has been decoupled using the 
BCS approximation 0 . The odd parity hybridization 
is such that, V(—r) = —V(r) in real space, or in fc-space 
V-k = —14 and mixes states with the same spin. 

The order parameters that characterize the supercon¬ 
ducting phase are, the inter-band superconducting order 
parameter given by, 

A a b = g a b ^ ^ (akcjb—k—a) ( 7 ) 

ka 

and the intra-band one, 

A bb — 9bb ^ bkgb—k—g )» (8) 

ka 
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where g a b and gbb are the attractive interactions. Al¬ 
though there is no attractive interaction in the a-band, 
we will investigate the existence of induced superconduc¬ 
tivity in this band. For this purpose we define the In¬ 
dependent anomalous correlation function in the a-band: 

A aa (fc,cr) = (dkaa-k-cr)- (9) 

This anomalous correlation function, as we will show be¬ 
low, turns out to be finite even in the absence of interac¬ 
tions in the a-band, due to the influence of hybridization 
and/or inter-band interactions. 

We use the equation of motion method to find the 
relevant Green’s functions and use the fluctuation- 
dissipation theorem to obtain from them the anomalous 
correlation functions above. Since we want to calculate 
the chemical potential, we also need to obtain the cor¬ 
relation functions (aj :a aka) and (frf ka bka) that yield the 
average number of particles in each band. 

IV. CALCULATIONS 


where, 


w ^4) 

A ab 

-14 

0 \ 

a:„ 

(« + £*) 

A * bb 

V- k 

-Vk 

A b b 

( w - 4) 

—A ab 

0 

V-k 

-a: 6 

(w + ^fc) / 


(i5) 

Notice that, from this system of equations we cannot 
compute all the correlation functions initially desired. 
We also need to solve another closed set of equations 
that is obtained when we calculate the equation of mo¬ 
tion for the anomalous Green’s function ((b- k - a \b k a))- 
This new set of equations is given by, 


( %2 = (( a ka\b ka )) \ 

V2 = «4 fc _>L» 


o o 

22 = ((M4ct» 


1 

\u 2 = {{a}_ k _ a \b\ a )) ) 


Vo/ 


In the next subsections, we obtain the relevant corre¬ 
lation functions and the energy of the excitations. 


In this section, we calculate the Green’s functions nec¬ 
essary to find the intra-band and inter-band order pa¬ 
rameters, as well as, the occupation numbers in the a 
and b bands. 

The equation of motion for the anomalous Green’s 
function ((b ] _ k _ a \a\ a )) is given by, 

( w + e -k ) ((^—fc—cr l a fccr)} + ^ab ({ a kcr l4er)) 

+ Aj>fc((4 fco .|4 fc(T }) + ^-fc«al fc _>L» = °- 

( 10 ) 

This generates two new Green’s functions for which we 
write the equations of motion, 

( w — 4) (( a fc<r|4cr)) 

+Aa6«&L fc _>L>> - V k ((b ka \a L» = 1 (11) 

and 

4 - 4) «M4<t» ~ A ab((a_ fc _j4j) 

+A bb (k,a)((bl k _ <J \a{ <J )) - V k {{a k<J \a\ a )) = 0. 

( 12 ) 

Finally, we find the last Green’s function namely, 
((alk-CT^io-)}! c l° ses the set of equations, 

(w + e“ fc ) «4 fc _>L» 

-A* ab ((b ka \al)) + Vl k ((bl k _ a \a{ a )). (13) 
It is helpful to write these four equations in matrix form, 


D. 


kal 


yi = «&L fc _>L» 

zi = «m4<t» 
V u i = «4 fc _>L» J 


\ 

(!) 

) 

loj 


(14) 


A. Excitation Energies 


The excitation energies of the system are given by the 
poles of the Green’s functions. These poles are obtained 
from the equation det(D) = 0 given by, 

w 4 - 2 A k ui 2 +B k = 0 (17) 


with 
Afc = 

Bk = 


,a 2 I ,62 


*66 


+1 A a b\ 2 + \v k y 


(44 — |t4| 2 + |A ah | 2 ) + 4| Vfc |“ | Aa,b | 2 -Fe^ 2 1A^^ | 2 

(18) 


where we have used the antisymmetric property of the 
hybridization, V- k = —14 and that the band energies 
are symmetric, i.e., e_ fc = e k . 

We assume without loss of generality that the order 
parameters A ab and A bb are real. Since the hybridization 
14 has to be purely imaginary to preserve time reversal 
symmetry, a term —2e k ^Re(A a bV k Abb) that appears in 
B k turns out to be identically zero and has not been 
written above. 

Notice that when B k = 0 we find zero energy solutions: 

oj 2 (w 2 - 2 A k ) = 0 (19) 

These zero energy modes are associated with topological 
transitions in the system 00 as we will see further 
on in the text. 

Finally, the excitation energies are given by, 
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V. SOLUTION OF THE PROBLEM 


The inter-band order parameter is defined by Eq. [7] 
From the closed set of equations, Eqs. [UjJ we can obtain 
the quantity y\ and using the fluctuation-dissipation the¬ 
orem we get the first gap equation, 


= — y 

4 n 


kcr 


(wi tanh(/?Wi/2) —W2tanh(/3w2/2)) 


1 

9ab 

(44 + \ V k\ 2 + A lb) ( tanh(/?wi/2) tanh(/3w 2 /2) \ 


UJt—UJZ 


UJf +U>2 \ 

At T = 0 K, this becomes 

1 1 V > 1 

'< ' 47T 


u 1 


Cl > 2 


9ab 


kcr 


0J1+UJ2 


c a 


i+- 


-\Vk\ 2 + A2 al 


UJ 10 J 2 


( 21 ) 


.( 22 ) 


For the intra-band gap term we get, 


— = -V 

g bb 47r ^ 


{uj\ tanh(^wi/2)o;2 tanh(/3w2/2)) 


^a2 


/ tanh(/3wi/2) tanh(/3w 2 /2) \ 


,2_, ,2 l 

J 1 W 2 V 


Wl 


U>2 




At T = 0 K, this is given by. 

1 _ 1 r 1 

/111 Air • ^ 


9b b 


kcr L 


W1 + W2 


,_a 2 


1 + 


W1W2 / 


(23) 


(24) 


Linear terms in 14 do not appear in these equations due 
to the assumptions that both A ab and A bb are real and 
that the system has time reversal symmetry in the ab¬ 
sence of superconductivity. 

The above equations involve the order parameters A ab , 
A bb and the chemical potential p through the dispersions 
of the bands. A full solution to the problem requires a 
self-consistent solution of a system of equations involving 
the three variables, A Q b, A bb and p. The equation for the 
chemical potential is obtained from the conservation of 
the total number of particles n, 


n = n a + n b 

™ = E + ( 6 LM) • (25) 

kcr 


The quantities (aj' 0 .afe cr ) and {b\ a bka) are obtained from 
the associated Green’s functions and 22 , in Eqs. (full) 
and (TTB1) . respectively). 

The total number of particle is given by, 


n = 


1 

(2^ 


47T 

T 


(4=> + & F,b ) I 


(26) 


where kp, a and kp, b are the respective wave-vectors for 
the a and b bands. We will assume that these bands are 
homotetic, i.e., 


(27) 


where ej, = k 2 /2m a - The quantity a = m a /m b is the ra¬ 
tio of the effective masses. Notice that, kF, b = kF, a /\fa. 
The Fermi energy is given by, Ep = k 2 F a ^/2m a ^) and 
the total number of particles can be expressed as, 


1 

47rfc F,a / 1 + a 3/2 \ 

(2tt) 3 

3 V a 3 / 2 ) 


(28) 


Using Eq.[ 
number as, 


we write the equation for the occupation 


1 = 


v 3/2 


37T 


2 kp „ V1 + <a 3 / 2 


E 

ka 


2- 


UI 1 + UJ2 


(4 + 4) (4-4 + A lb - \Vk\ 2 ) + e-k A bb 


W1W2 (wi + OJ2) 


, (29) 


at T = 0 K. 

Eqs. HU HU and HU define our problem with intra and 
inter-band superconductivity. For solving them summa¬ 
tion in k is changed to integration using, 

< 30 > 

Here and below tilde quantities mean, for wave-vectors 
normalization by kF,a and for energies normalization by 
Ep■ V is the volume element that we take, without loss 
of generality, as V = 1. 

We will not investigate in this work the full self- 
consistent problem involving the two superconducting or¬ 
der parameters and the chemical potential. Instead, we 
consider the particular cases where we have, either intra 
or inter-band interactions only. Furthermore, we want to 
be able to extend our calculations to the strong coupling 
regime where g ab and g bb are very large. As usual, in 
the study of the crossover from the weak to the strong 
coupling reg imes, we introduce two scattering lengths a s 
and a sb [32|, |33( which replace the coupling constants g ab 
and g bb , respectively. Since the integrals are done for all 
values of k , they diverge and have to be regularized. We 
employ here the standard procedure to eliminate these 
ultra-violet divergences (see below). 

• Inter-band case 

In the pure inter-band case g bb is zero and there are 
only inter-band attractive interactions in the sys¬ 
tem. We are left with two equations for A ab and the 
chemical potential p to be solved self-consistently. 
These are given by, 


1 + a 

47T 3 



1 

Cb\ + UJ 2 


2 

(1 + a)k 2 


(a 2 6 + 44 + I^1 2 ) 

W1W2 (d>i + W 2 ) 


1 


kpa s ’ 


Cfc = efc - p 

4 = ae k - Mi 


(31) 
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and the number equation, 


VI. NUMERICAL SOLUTION 


A/2 


167T 2 \ 1 + a 3 / 2 

(3+3) ( 


d 3 k 


2 - 


~a , 

£ fc+ e fc 

Wi + W 2 


£ fc £ H^| 2 - 


-|A 


a6| 


W 1 W 2 (wi + w 2 ) 


= 1 . 


(32) 


We used the definition of the scattering length for 
the low energy limit of the two-body problem in the 
vacuum [32j . 

— = - -b [ d 3 k [/1 (k) - fi (k -a 00)] (33) 

Qab 47 TCL S J 


A. Inter-band case — gbb = 0 


In this section we solve self-consistently Eas.l3Tlandl5!21 
for the inter-band superconducting order parameter A ab 
and the chemical potential /i. We consider a three dimen¬ 
sional system and take for the hybridization the form 
V(k) = i'y (k x + k y + (3k z ). For /? < 1 we want to de¬ 
scribe a tetragonal system where the hybridization is 
smaller between planes. Since we are here interested in 
studying the effect of hybridization on superconductiv¬ 
ity, we take the quantity l/kFd s negative and small, such 
that the system is in the weak coupling BCS regime. The 
results for A ab and g are shown in Fig. ©■ They corre¬ 
spond to fixed l/kpo-s = —0.5 and a ratio for the masses 
a = 0.1. Also the parameter /? = 0.1 in 14. We see 


where a s is the s-wave scattering length, m* = 
a/(\ + a)m a and a = m a /mb is the ratio of the 
effective masses of the a and b quasi-particles. The 
function f\{k) in this case is given by, 

f 1 (k)= 1 | ( e fc4 + l+fc| 2 +l A a&| 2 ) 

UJ!+UJ2 WiW 2 (Wi + W 2 ) 


Since the integral extends to infinity the subtrac¬ 
tion of the last term on the right hand side of Eq. 
(EH) regularizes the ultra-violet divergence in this 
expression. 

• Intra-band case 

In this case, we take g a b = 0 and the gap equation 
is given by, 


1 

k F a sb 



1 

Wl+W 2 


U)lU>2 


1 

ah 2 


(34) 


where we replaced gbb by the scattering length. The 
occupation number equation is now given by, 


3 / a 3 / 2 \ 

16-7T 2 v 1 + a 3 / 2 ) 


d 3 k 



c k~ 


+ U)2 


e a k \A bb \ 2 (3 + 3) (33-|Vfc| 2 ) 

0J10J2 (wi+w 2 ) W1W 2 (wi + w 2 ) 


1.(35) 
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(b) 


Notice that in each of the cases above the excitation 
energies u)i and w 2 are different. They are obtained from 
Eqs. [20] making either A ab = 0 or A bb = 0, as appro¬ 
priate. As before tilde quantities mean, for energies nor¬ 
malization by the Fermi energy Ep, for wave-vectors by 

kF,a- 


FIG. 1. (Color online) The normalized chemical potential 
ft and the inter-band superconducting order parameter A ab 
as functions of the strength of the hybridization for fixed 
1 /kFd a = —0.5. This value for the coupling is typical of the 
weak coupling BCS regime. We used a = 0.1 for the ratio of 
the effective masses and /3 = 0.1 for the anisotropy parameter. 
Notice that the hybridization induces a BCS-BEC crossover. 
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from Fig. m that the normalized chemical potential jl 
decreases as the strength of the antisymmetric hybridiza¬ 
tion 7 = 7 /vf increases. This decrease is accompanied 
by an increase of the inter-band superconducting order 
parameter A a b- There are at least three very interesting 
points to be noticed in these figures. First, the enhance¬ 
ment of superconductivity by antisymmetric hybridiza¬ 
tion which had already been noticed [0|. Second, the 
combined behavior of the chemical potential decreasing 
7 and the concomitant increase in A a b is a clear signa¬ 
ture of a BCS-BEC crossover induced in this case by an 
increase in hybridization [l3j . Finally, notice the discon¬ 
tinuous behavior of ft and A a b for 7 « 0.65, due to a 
topological phase transition associated with the appear¬ 
ance of gapless modes in the spectrum of excitations for 
this value of 7. This is due to a vanishing of the quan¬ 
tity Bk in Eq. m • For Abb = 0, the condition Bk = 0 
implies 14 = 0 and e k e b k + A 2 b = 0. The intersection of 
these two surfaces give the line of zero energy excitations 
associated with the topological transition. This transi¬ 
tion is very sensitive to the ratio of the effective masses 
and for a = 0.5 it has disappeared for the range of 7 in 
Fig. m 


B. Intra-band case — g a b = 0 


Here we discuss the pure intra-band case, where g a b = 
0 in Eq. 0 In this case superconductivity is character¬ 
ized by the intra-band order parameter Abb, defined in 
Eq. [HI Taking A a b = 0 we are left with a system of two 
equations to be solved self-consistently. As in the pre¬ 
vious inter-band case, we assume for the antisymmetric 
hybridization the form V(k) = i'y ( k x + k y + /3k z ) where 
14 is purely imaginary to preserve time reversal symme¬ 
try. The self-consistent solutions for the order parameter 
Abb and the chemical potential in terms of 7 /vf =7 are 
shown in Fig. [2] of this section. Here we recover the usual 
deleterious behavior of the hybridization on intra-band 
superconductivity Si- We can see that Abb differently 
from A ab decreases with increasing 7 eventually vanish¬ 
ing at a superconductor quantum critical point. 


VII. INDUCED SUPERCONDUCTIVITY 


An interesting result of our investigation is that, in 
spite of the absence of interactions in the a-band, the 
anomalous correlation function A aa (k) may assume fi¬ 
nite values, signaling the presence of superconducting 
correlations in this non-interacting band. The anomalous 
correlation function A aa {k) = (ako-a-k-cr) can be easily 
calculated from the associated Green’s function (iq) and 
we get, 



FIG. 2. (Color online) The normalized chemical potential 
ft and the intra-band superconducting order parameter Abb 
as functions of the strength of the hybridization for fixed 
1/kFdsb = —0.5, typical of the weak coupling BCS regime. 
We used a = 0.1 for the ratio of the effective masses and 
P = 0.1 for the anisotropy parameter. 


A aa {k) = - 


1 

4-7T 


2e b k VkA a b Abb (1141 2 A 2 ab ) 


x 


/ tanh(/3wi/2) tanh(/ 3 w 2 / 2 ) \ 




w 1 


W 2 




(36) 


At T = 0 K, this becomes, 


- 4?r 


2e b k VkA ab - Abb QVk[ 2 - A 2 ab ) 

W1W2 (UJl + 0J2) 


.(37) 


This induced anomalous correlation function in the a- 
band arises due to the influence of hybridization and/or 
inter-band interactions. A linear term in 14 now appears, 
indicating that A aa {k) can be complex. This correlation 
function is a sum of antisymmetric and symmetric con¬ 
tributions. Notice that e k , w 1 and W 2 are even functions 
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of k. The antisymmetric term proportional to 14 corre¬ 
sponds to induced p-wave anomalous correlations due to 
the odd-parity of 14- It preserves time reversal invari¬ 
ance since A ab is real and 14 is purely imaginary. This 
term corresponds to the component in = 0 of the l = 1 
p -wave state. 

Eq. (1371) for the induced anomalous correlation func¬ 
tion A aa is very interesting, and deserves special atten¬ 
tion. We notice that there are three possible situations 
involving the relative magnitudes of the order parame¬ 
ters, that will be discussed now. 

1. The first situation is met when the second term 
in the right side of the above equation can be ne¬ 
glected, i.e., when 2e^|I4|A ab > A bb (\V k \ 2 - A 2 ab ) 
with the ^-dependent quantities calculated at kp. 
Observe that this is naturally true when A bb is very 
small (compared to A ab ) and also when |T4| ~ A ab , 
regardless of the value of A bb . Then we can write, 


A aa (k) 


1 

2n 


e b k V k A ab 

0J\U}2 (Wi + OJ2) 


(38) 


Since the parity of A aa follows that of V k which is 
antisymmetric, A aa (k) has odd-parity in k- space. 
Then the superconducting fluctuations induced in 
the a-band have a p-wave character. 

2. The second one happens in the opposite limit, that 
is A bb (| T4 1 2 - A 2 ab ) > 2e||14|A a6 , which can be 
obtained, for instance, for a vanishingly small V k , 
resulting in, 


A aa (fc) 


1 

47T 


A bbA 2 ab 

OJ1UJ2 (wi + w 2 ) 


(39) 


See that in this case A aa is completely symmetric. 

3. The third case is obtained when A ab ss 0, 


■<*> = s 


—A bb \V k \ 2 

OJ 1 UJ 2 (wi + w 2 ) 


(40) 


This corresponds to the more conventional induced 
superconductivity, as in the proximity effect. It is a 
second order effect in the hybridization. In the last 
two cases the induced superconductivity is clearly 
s-wave. 


VIII. MAJORANA FERMIONS 

There is nowadays a great interest in the study of Ma- 
jorana fermions as possible candidates for use as q-bits 
in quantum computers. In an exciting paper jl5| . Kitaev 
has shown that one-dimensional p-wave superconductors 
can exhibit a non-trivial topological phase with Majorana 
fermions at the ends of the chain. A great effort is being 


done to implement this idea of Kitaev in an actual physi¬ 
cal wire [341 ]. The main difficulty is of course that p-wave 
superconductors are scarce in nature and to circumvent 
this, many suggestions have been proposed. The most 
successful one combine a mixture of spin-orbit interaction 
and external magnetic field to induce in a semiconductor 
wire p-wave type of superconductivity [34j. 

The results above suggest a new mechanism for induc¬ 
ing p-wave superconductivity in a wire and consequently 
to obtain Majorana fermions. The idea is to consider the 
Hamiltonian, Eq. [6l as an effective Hamiltonian for the 
following system. A bulk BCS superconductor with a b- 
band of electrons with an attractive interaction responsi¬ 
ble for the superconductivity characterized by the s-wave 
order parameter A bb . On top of this material, there is 
a normal wire with a non-interacting band of a-electrons 
which are coupled to the bulk substrate through an an¬ 
tisymmetric hybridization and a weak attractive interac¬ 
tion g ab . The latter gives rise to the inter-band pairing 
A ob among the electrons of the wire and the substrate. 
Then, under the conditions 1 of Section VI, we obtain 
an induced p-wave superconductivity in the wire. 

In order to verify if the induced superconductivity in 
the wire gives rise Majorana fermions at its ends, we 
need to write the wire Hamiltonian in terms of Majorana 
operators. First, we present the Hamiltonian of the wire 
in terms of fermions operators, 


‘ ‘ ii'/ rr — 


N 

Y (^t u ) u 3 + H d )dj^ 

j =1 

-* Y • v 1 "/H• 1 • 

3= 1 

N -1 

+ Y A ( U i d 3 + 1 “ U 3+l d j + d )+l U ) - d ) U )+l) >( 41 ) 


N-l 


i =1 


where u' and dl are creation operators for spins up and 
down in the a-band of the wire, respectively. The spin de¬ 
pendent chemical potentials are given by, /x-m) = H ± h, 
where h an external magnetic field. The quantity t is 
a nearest neighbor hopping and A = | Ay | is the in¬ 
duced antisymmetric pairing. The antisymmetric na¬ 
ture of this coupling has already been taken into account 
when writing the Hamiltonian in the form above. Notice 
that A ij = ^2 k A aa {k)e lkij ' i ~ rj which is antisymmetric 
in real space since A aa (k) oc 14 is the induced p-wave 
superconductivity in the wire due to the antisymmetric 
mixing between the electrons in the wire and the bulk 
( V- k = —14). This pairing interaction couples electrons 
of opposite spins. It corresponds to the mi = 0 compo¬ 
nent of the l = 1 p-wave state. This differently from the 
m = ±1 components does not break time reversal sym¬ 
metry. This is broken in the Hamiltonian above by the 
external magnetic field h. 

We can rewrite the fermion creation and annihilation 
operators of Hamiltonian Eq. |TT] in terms of Majorana 
fermions operators. These are given by the following def- 











initions 


? 


Uj — (CXBj + iCXAj ) 

dj = —^= {pBj + ipAj) 


(42) 


and similar equations for their complex conjugates, notic¬ 
ing that the Majorana particles are their own antipar¬ 
ticles, i.e., a' A - = ctAj, a' gj = a Bj, = Paj and 

(3] jj = (3 b j- In the Majorana basis the Hamiltonian 
Eq. [41] becomes, 


N 

Hwire — 7^ ^ ^ (/R 'H-Bj^Aj T P[ PBj PAj ) 

3 =i 


N -1 

— * E t (aBjO-Aj+1 — OtAjOCBj+l + pBjPAj+1 — PAjPBj+l) 
3 =1 


AT-1 

+* E A (otAj(3Bj+l+CXBj/3Aj+l — ClA}\-lpBj — ^Bj+lPAj) ■ 

3=1 


(43) 


We introduce new hybrid Majorana operators, like in 
Ref. [Uj to get, 

7a j = a Aj ± Paj 

1Bj=aiBj±pBj, (44) 

and rewrite Eq. 03] in term of these new operators. We 
obtain, 


(45) 

This equation shows that the magnetic field h couples the 
operators 7 . When this field is zero, the Hamiltonian, 
Eq. 05] reduces to that of two decoupled Kitaev chains. 
It is easy to see that for h = 0, this model has a topolog¬ 
ical superconducting phase for \p/2t\ < 1 as the Kitaev 
model |l5|. Let us consider p = h = 0 and A = t, such 
that the system is in the topological phase. In this case 
the Hamiltonian of the wire is given by, 

N 

H wire = -it E (lBj7A3 +1 - 7Aj7Bj+l) ■ ( 46 ) 

3=1 

Notice that the hybrid Majorana operators and 
do not enter in this Hamiltonian, so that, there are two 


it _ 


Ifl 

T 


N r 


E ( 7 ^ 7 % +7 Bj 7 Aj )- y (lijlAi+lBjlX 


3=1 


N 


i(t — A) y- / 

E( 

3=1 
N 

zt 


2 

i(t + A) 


7bj 7 ) 47+1 7Aj7Bj+i 


) 


7Bj7Aj+l 7Aj7Bj+l 


3=1 


unpaired hybrid Majorana fermions in the left end of the 
chain. The same occurs for 7 AN and j gN on the right 
end of the chain. We can combine these four Majorana 
fermions as two ordinary fermions, ipi = ( 7 X 1 + *7ai) / 2 
and ipN = (7bn ~ * 7 an) / 2 - Also we can combine the 
Majoranas at different ends of the wire, such that, 3p + = 

(7m + Han) / 2 and = (7bn ~ *7ai) / 2 - Fi g- El 
shows schematically these combinations. Next, we de- 

v + 



Mb 

FIG. 3. (Color online) Schematic figure showing the two Ki¬ 
taev chains, + and —, associated with the wire. There are 
four Majorana states at the ends of the wire. The local and 
non-local combinations of the Majoranas in fermion states are 
also shown. 

fine local pseudo-spin operators as in Ref. [25j . 

S 1 (JV) = -j 7 l(N)i a i 37 l{N)j, ( 47 ) 

where the components of <7 are the Pauli matrices, i 
and j can be + and —, and 71+ = 7/y 71 _ = 7 A1 , 
7 n+ = 7an and 7iv- = 7bn- Using the anticommuta- 
tion relations of the Majorana operators, {7i+,7,_} = 0 
and {7 i+ (_), 7j + (_)} = 2 5ij, we obtain s \(n)- It turns out 
that both si and s n have only the y component different 
from zero. In this sense we conclude that these composite 
Majorana fermions behave as Ising spins. Furthermore, 
we have, spN)y = (*7i(at)-7i(jv)+)/2, that can be rewrit¬ 
ten as, 

s l(N)y — ~ 1W)i (48) 

where 'ipum were defined above. It is interesting to write 
this expression in this way because we can relate these op¬ 
erators to the occupation numbers of the fermion states 
at each end of the wire. 

As in Refs. [2fJ, we connect our system (supercon¬ 
ductor bulk and wire) to a capacitor as in Fig. 0] to 
control the occupation numbers of the fermion states 
in the ends of the wire. There will be a total of 
four occupancy states: |0)i|0)at, |1}i|1)at, |0)i|l)jv 

and |1)i|0)jv- Notice that the first two states have 
an even number of fermions and the last two an 
odd number. We can define two density matrices: 
Podd = 1/2 (|0) r 11 )jv i (0 |jv( 11 + 11)i|0 )jv i (11 iv (01) 

with odd parity states and Peven = 
1/2 (|0)i|0)jv i(0|tv(0| + |1) 1 |1)jv i<1|jv<1|) with even 
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parity states. We then calculate the value of the 
correlation function < Si y S;vy > in these different parity 
states. Since, < A >=Tr pA, we get, 

('N y^Ny) odd = {^lySiN) even — (49) 

This shows that there are long range correlations between 
the end states of the wire, which are independent of its 
length. Controlling the charge of the end states we can 
vary the sign of these long range correlations. These 
highly non-trivial properties of this type of p -wave su¬ 
perconductor in its non-trivial topological phase can cer¬ 
tainly find useful applications, for example, in quantum 
computers. 



FIG. 4. (Color online) Schematic figure showing the bulk- 
wire system connected to a capacitor to modify the fermion 
number via a charging energy: U(n) = Q 2 /(2C) — QV(t). 
Here, the number of fermions is changed by variating V I2~J |. 


IX. CONCLUSIONS 

In this paper, we have studied the influence of an odd- 
parity hybridization in the superconducting properties of 
a two-band system with intra (in only one of the bands) 
and inter-band attractive interactions. First, we dis¬ 
cussed the origin and conditions for antisymmetric hy¬ 
bridization, and showed that it occurs in many systems 
of great interest as organic materials, copper oxides and 
heavy fermion systems. In these systems the relevant 
bands consist of electrons whose angular momenta dif¬ 
fer by an odd number, as for s-p , p-d, and ((-/-electrons, 
which is the required condition for antisymmetric mixing. 
Antisymmetric hybridization is also a crucial element for 
the appearance of topological insulating phases in multi¬ 
band systems [30| . 

Using a Green’s function approach, we obtained a sys¬ 
tem of three self-consistent equations for the intra and 
inter-band superconducting order parameters and the 
chemical potential. Our approach incorporates a scat¬ 
tering length and allows to treat the cases of weak and 
strong inter and intra-band attractive interactions and 
consequently to study the crossover between the weak 
coupling BCS regime and the Bose-Einstein condensa¬ 
tion (BEC) of pairs. This crossover has been intensively 
studied theoretically and its experimental observation 


was only possible with the advancement of cooling tech¬ 
niques for atomic gases. In these systems it is possible 
to control the strength of the many-body interactions, 
which is not feasible in condensed matter systems. Here 
we obtained the remarkable result that hybridization can 
drive the BCS-BEC crossover even keeping the interac¬ 
tions fixed at values characteristic of the weak coupling 
regime. Since hybridization can be tuned in solid state 
matter by applying pressure, such as to vary the overlap 
of the wave-functions, this opens a unique possibility for 
investigating experimentally the BCS-BEC crossover in 
condensed matter systems. 

We have studied the induced superconductivity that 
appears in the non-interacting band of our two-band sys¬ 
tem due to the combined effects of hybridization and 
inter-band interactions. We have shown that in the case 
the former is antisymmetric, the induced superconduc¬ 
tivity has a p -wave character. This induced pairing cor¬ 
responds to the to; = 0 component of the l = 1 p -wave 
state and does not break time-reversal invariance. 

Our two-band model represents an effective model for 
a system consisting of a bulk superconductor, described 
by the 6-band, and a normal wire deposited on top of 
it with a non-interacting band of electrons (the a-band). 
The normal wire is coupled to the bulk substrate by an 
antisymmetric hybridization and an attractive inter-band 
interaction. We have shown that one component of the 
induced superconductivity in the wire due to hybridiza¬ 
tion is of the p -wave type, more specifically it is associ¬ 
ated with the time reversal invariant projection, to; = 0 . 
We have studied the excitations in a wire with this type 
of induced superconductivity using a Majorana represen¬ 
tation. We found that it has a topological superconduct¬ 
ing phase that supports four Majorana modes, two at 
each end of a finite chain. This is different from Kitaev’s 
spinless chain model that corresponds to the to; = ±1 
p -wave state and has a single Majorana at each end of 
the chain in its non-trivial topological superconducting 
phase. Although two Majoranas give rise to a conven¬ 
tional fermion, we have shown using results of previous 
works that the charge states of the fermion modes at the 
ends of the to; = 0 , p -wave superconducting chain have 
highly non-trivial long range correlations. Since these 
states can be controlled this system may have interesting 
properties for applications. 
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